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Abstract 

We discuss some aspects of a new noncombinatorial fermionic approach 
to the two-dimensional dimer problem in statistical mechanics based on the 
integration over anticommuting Grassmann variables and factorization ideas 
for dimer density matrix. The dimer partition function can be expressed as 
a Gaussian fermionic integral. For regular lattices, the analytic solution then 
follows by passing to the momentum space for fermions. 

1 Introduction 

The dimer problem was originally formulated to describe the entropy of diatomic 
molecules adsorbed on a crystal surface [1]. In a more abstract sense, the two- 
dimensional (2D) dimer model is one of only few nontrivial exactly solvable models 
in statistical mechanics [2-8]. In its original state, this is a purely combinatorial 
problem, close in its inherent spirit to the 2D Ising model [4] . The exact solution for 
the closed-packed 2D dimer problem on the standard rectangular lattice was first 
obtained by a direct combinatorial method [2,3]. A number of more complicated 
dimer lattices have also been analyzed [5,6]. Traditional combinatorial approaches 
of this kind, however, are somewhat complicated, and are different, in their spirit, 
from the methods commonly used in condensed matter physics. In this report we 
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comment on a simple fermionic approach to the 2D dimer problem based on the 
integration over anticommuting Grassmann variables (nonquantum fermionic fields) 
and the mirror-ordered factorization procedure for the dimer density matrix [7,8]. 
Neither transfer-matrices nor combinatorics are used. The factorization method 
rather resembles the idea of insertion of the Dirac unity, E|o)(a| = 1, when mak- 
ing transformations in quantum mechanics. We first reformulate the combinatorial 
dimer partition function, Q, in terms of the commuting nil-potent (^-variables, 
and then introduce into the resulting expression the anticommuting (a)-variables 
by factorization of the local bond Boltzmann weights. Eliminating the ^-variables 
in the resulting (rj, ^-representation, we come to a purely fermionic expression for 
the same partition function Q. Even for the inhomogeneous distribution of the 
dimer weights over the lattice, the partition function is expressible as a Gaussian 
fermionic integral. Equivalently, the closed-packed 2D dimer problem is a theory of 
free fermions on a lattice. For regular lattices, the analytic solution then follows by 
transformation to the momentum space for fermions. 



2 Grassmann variables 

Let us remember first the basic rules of anticommuting analysis [9]. Grassmann 
variables are classic (nonquantum) fermionic numbers purely anticommuting to zero. 
Given a set of Grassmann variables a±, a^, we have: aidj + ajdi = 0, = 0. The 
first important identity is that for the product of linearly transformed variables: 

N 

bi b 2 .... b N = det A ■ a\ a 2 ...a N , h = ^2 ^ij a J > (^-1) 

3=1 

where A is the matrix of transformation a — > b. The determinant here appears due to 
the known interrelations between fermionic algebra and determinant combinatorics. 
The rules of fermionic integration were first formulated by Berezin [9]. For one 
variable, we have: 

J da j ■ dj = 1 , J da j ■ 1 = . (2.2) 

In the multiple fermionic integral the differential symbols dai, da^ are again 
anticommuting with each other and with the variables [9]. Due to the property 
a* = 0, any natural function defined on the set a±, . . . ,a,N can be written as a finite 
polynomial in these variables: 

/ (ai, a 2 , a N ) = f + fj a i + ' ' ' + /123...JV aia 2 ...a N , (2.3) 

(i<i<A0 
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where f , . . . , fu3...N are numerical parameters. Integrating a polynomial function 
like (2.3) according to (2.2), we obtain: 

J da N . . . da 2 da 1 f (a^ a 2 , ■ ■ ■ , a N ) = f 12 3...N , (2.4) 

where fm...N is just the coefficient in the last term of the polynomial (2.3). We see 
that fermionic integration is a trivial problem if the integrand function is already 
known in a polynomial form. However, this may be not the case in applications. 
From (2.1) and (2.4), one can deduce the rules of change of variables in a fermionic 
integral under a linear substitution. As compared with the rules of a proper commut- 
ing analysis the only difference is that the Jacobian will now appear in the inverse 
power. 

An important role in applications play Gaussian fermionic integrals [9,10]. The 
Gaussian integral of the first kind is related to the determinant: 




where {aj, a*} is a set of purely anticommuting Grassmann variables, the ma- 
trix A is arbitrary. By convention, the variables and a* can be considered as 
complex-conjugated fermionic fields, otherwise these are independent variables. The 
fermionic exponential is assumed in the sense of its series expansion, the series termi- 
nates at some stage due to the property aj 2 = 0. The correspondent finite polynomial 
also can be obtained by multiplying elementary factors exp (diA^a*) = 1 + diA^a*. 
The Gaussian integral of the second kind, for real fermionic fields, is related to the 
Pfaffian: 

f (i n n \ 

I daj^ ... da 2 da\ exp - S ai Ai o a i = Pfaff A > A + A T = 0, (2.6) 

J v 2 i=ij=i / 

where the matrix A is now skew-symmetric (Aij + Aji = 0, An = ). The Pfaffian 
is some combinatorial polynomial in elements A^ well known in mathematics for 
a long time. In fact, the Pfaffian combinatorics is identical with that of fermionic 
Wick's theorem. The Pfaffian and determinant of the associated skew-symmetric 
matrix are algebraically related: det A = (PfaffA) 2 . This identity can be most 
readily proved just in terms of integrals like (2.5) and (2.6), making use of a suitable 
change of variables in (2.5). 
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3 The dimer problem 



Let us consider the closed-packed 2D dimer model on a rectangular lattice net with 
the most general inhomogeneous distribution of dimer weights [7,8]. The dimers 
are objects living on lattice bonds. A given bond may be either free or covered by 
a dimer together with the two adjacent lattice sites. The closed-packing condition 
means that each lattice site must be occupied by one and only one dimer. The 
lattice must have an even number of sites to be covered completely by dimers in 
a closed-packed fashion. An example of a closed-packed dimer configuration on a 
rectangular net can be seen in Fig. la in [7]. 

Let the lattice sites be numbered by the integer Cartesian coordinates mn, where 
m — 1, . . . , M, n = 1, . . . , N run in horizontal and vertical directions, respectively. 
MN is the total number of sites. We define t$ n and t$ n to be the dimer weights 
(activities) for horizontal (mn \ m+ln) and vertical (mn | mn+1) bonds, respectively. 
The weight of a free bond is 1. The dimer partition function is: 

M N (i) ( 2 ) 

Q=E(nnfi 4 }. (3- 1 ) 

{jV} m=ln=l 

where and JV^ are dimer occupation numbers taking the values 1,0, depending 
of whether a given bond is covered or free, and the sum is taken over all closed- 
packed configurations. Notice that in a completely homogeneous and symmetric 
case, t$ n = t$ n = 1, the partition function yields simply the total number of closed- 
packed dimer configurations for a given lattice. The thermodynamic interpretation 
is only possible if there are few sorts of dimers with different weights (energies). 

We need to transform Q into a fermionic Gaussian integral. The combinatorial 
definition (3.1) is not very suitable in this respect since the symbol of the sum in (3.1) 
is only a conventional prescription for counting over closed-packed configurations. 
In fact, we need some kind of a formal averaging which will select the closed-packed 
configurations automatically. This can be realized in terms of the commuting nil- 
potent variables [7]. With each lattice site we now associate the commuting nilpotent 
variable i] mn , such that n^ n = 0, and write: 

M N 

Q = S P { II n^+tmnVmnVm+l n) (1 ^trnnVmnVmn+l ) } , (3-2) 
(rj) m=ln=l 

with the averaging rules for one variable: 

SP (l\Vmn\V 2 mn\V 3 mn\-)= ( | 1 | | | ...) ■ (3.3) 
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The global ^-averaging in (3.2) assumes local averagings like (3.3) at all sites. In (3.2) 
we also assume the free-boundary conditions: % +ln = i] m ^ + i = 0. Expression (3.2) 
exactly reproduces the combinatorial partition function (3.1). The product of the 
Boltzmann factors 1 +ti]i]' forming the density matrix in (3.2) generates all possible 
dimer coverings of a lattice by 'dimer molecules', rjtr)', with no restrictions, while 
the 77-averaging according to (3.3) selects the closed-packed configurations. These 
are the configurations in which each lattice site is covered by one and only one 'atom' 
or 'monomer', r\ mn , which can only be a part of some dimer. 



4 Fermionization 

Let it be given two Grassmann variables, a and a*. The elementary Gaussian expo- 
nential is e Xaa * = 1 + Xaa*, where A is a parameter. Noting (2.2), for the associated 
Gaussian averages we then find: / da*dae Xaa * {1, a, a*, aa*} = {A, 0, 0, 1}. Making 
use of these rules, we can factorize the weight factors from (3.2) as is shown in (4.1) 
below. For the whole lattice, we introduce a set of purely anticommuting Grassmann 
variables, { a mn , a* mn , b mn , b* mn }, and write: 

1 trnnVmnVm+ln 

= I da* mn da mn e a ™<n(l + t« a mn r] mn ) (1 + a* mn r} m+ln ) ^ A1(£) 
= Sp { A mn A* m+ln } , 

y&mn ) 
1 trnnVmnVmn+1 

= I db* mn db mn e & ™« & ™(l + t$ n b mn r] mn ) (1 + b* mn r] mn+l ) ^ ^ 

= Sp { B mn B* mn+l } . 

\Pmn) 

In the last lines we introduce the abbreviated notation for the arising factors (to 
be called shortly Grassmann factors), A mn , B mn , A^ +ln , B^ n+1 , while Sp (...) stands 
for the symbol of the Gaussian fermionic averaging. Neglecting the averaging sym- 
bol, the Boltzmann weights from (3.2) are now presented in a factorized form: 
A mn A* m+ln , B mn B^ n+1 . In general, for the whole lattice, there are four Grassmann 
factors, A mn , A* mn) B mn , B^ n , which all involve the same variable i] mn associated 
with a given mn site. The idea of the next step is to place nearby the above four 
factors and to average over r\ mn in each group independently thus passing to a purely 
fermionic expression for Q. The obstacle to this method is that the individual Grass- 
mann factors are neither commuting nor anticommuting with each other. It might, 
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therefore, be difficult, in general, to find the four relevant factors nearby in a global 
product. The problem of a suitable ordering of the non-commuting Grassmann fac- 
tors thus arises. This ordering problem can be solved applying the mirror-ordered 
factorized procedure [7]. The related ideas were first developed in the context of 
the 2D Ising model [11,12]. The result is that the density matrix from (3.2) can be 
represented in the following 'mirror-ordered' factorized form [7,8]: 



N M , M 



q(v) = s P { n [ n B ^n ■ n ^An^]}. (4.2) 

(a, O) n=l m=l m=l 

The 77-averaging can be performed at the junction of the m-ordered products in (4.2) 
[7]. This results in a purely fermionic expression for the partition function, which 
can in turn be simplified by integrating out the extra fermionic variables. In this 
way we obtain [7]: 

. N M ™_ M N 

Q = J UU dc mn exp{ ]T ]T [tVc mn c m+ln + (-l) m+ HWc mn c mn+1 ] } , (4.3) 

n=l m=l m=l n=l 

with CM+in — CmN+i — 0, where c mn are new totally anticommuting Grassmann 
variables. The product in the measure is ordered as follows. First, we multiply dc mn 
over m, with fixed n, then we multiply the resulting products over n. In the field- 
theoretical language, the fermionic form in the exponential is called action, since the 
action is quadratic in fermions, we deal here with a free-fermion field theory on a 
lattice. It is also possible to eliminate the sign factor (— l) m+1 from the action. This 
can be done by the rescaling of the variables in the integral: c mn — > c mn i m +3//2 , 
dc mn — > dc mn i~ m ~ 3//2 . This yields equivalent representation [8]: 

N M ^ M N 

Q = J II II dc ™ exp { ^2Y1 [tmnCmnCni+ln + % ^C mn C mn+ i ] } . (4.4) 
n=l m=l m=l n=l 

The dimer partition function is now expressed as a simple fermionic Gaussian inte- 
gral. The exact solution for a regular dimer lattice can now be obtained by passing 
to the momentum space for fermions (Fourier substitution). This is illustrated in 
the next section with examples of regular (homogeneous) rectangular and brick- 
hexagonal dimer lattices. It is interesting that these two models exhibit quite dif- 
ferent properties. While there is no phase transition for the standard rectangular 
dimer lattice, there is an exotic phase transition with frozen low-temperature phase 
for the hexagonal lattice. 
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5 Analytic results for regular lattices 



Let us consider the 2D dimer model on the standard homogeneous rectangular lat- 
tice. Assuming t$ n = t±, t$ n = t 2 in (4.4), the partition function is: 



N M ™ M N 



Q = J II II dc mn ex P E E [tlC mn C m+ln + it 2 C mn C mn+1 ] , (5.1) 
n=l m=l m=l n=l 

where c^+in = 0, c m ^+i = 0. We deal here with a finite lattice with a free boundary 
of size M x N, with MN sites. In what follows we assume that M and N both are 
even. It is interesting that integral (5.1) can be performed exactly despite of that 
the translational invariance is broken by free boundary. The diagonalizing Fourier 
substitution is [7,8]: 



C-m.n. 



2i m+n+ l ™" ( npm \ . / nqn 



^(M+lXJV + l) KM + l > KN + l 

In the momentum space, we find [8]: 

N M M N 

EE 

p=l q=l 



f N M ^ M N np nq 

J o=lp=l p=lo=l M+1 iV i- 1 



where Cpg, c P9 are the new variables of integration, p = M+1— p, q = N+l—q 
are 'conjugated' momenta selected by the orthogonality relations for the Fourier 
eigenf unctions from (5.2). Notice that the variables c P g, c pq enter, simultaneously, 
into the pq and pq terms of the total sum in (5.3) (in the given case, these terms 
appear to be equal to each other). To single out explicitly the true independent 
variables (which is needed to perform the integral) we have to combine these pq and 
pq terms together and then to reduce the sum to a half-interval with respect to M or 
N (choice between M or N is here arbitrary) in order to avoid the double counting 
of the same terms. Making the reduction with respect to M, the partition function 
becomes: 

r iM N w N ^ 

Q = / II II dc pq dc PQ ex P E E [ c m c pi ( 2t i cos 77—7 + 2it2 cos tttt) \ ■ ( 5A "> 

J p=lq=l p=lq=l JW + 1 JV "T 1 

The value of the integral readily follows from elementary rules like (2.2). In the 
resulting product we then make the q <-> q (or % <-> —%) symmetrization, and find: 



Atj cos 2 * P + 4i| cos 2 ' f/ 



Q=nn 

p=l q=l 



M+1 z N + l 



evenM,iV. (5.5) 
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This is the exact solution for the 2D dimer model on a finite rectangular lattice with 
free boundary. This result was first obtained by the combinatorial method [2,3]. 

Taking the limit of infinite lattice, we obtain the free energy per site for the 
standard rectangular dimer lattice: 



1 



1 



2 2 



-/3/ REC = ^lnQ mn ^ = -^J J dpd q ln[4tl C os 2 p + Ulcos 2 q\ . (5.6) 



The free energy (5.6) can be represented in other equivalent forms, in particular, it 
can be written in the form [8]: 

-/^ilnfMj + f + i//^, (5.7) 

*2 *2 

where G = 0.915 965 594 is Catalan's constant. The term G/n corresponds to a 
purely combinatorial contribution to the free energy (and entropy) while the ther- 
modynamics is governed by the doubled-integral term. In thermodynamic interpre- 
tation, the dimer weights are to be chosen in the form t a = exp (e a ), £ a = (3E a , 
a = 1,2, where E a are the dimer energies (more precisely, the true energies are —E a ), 
and P = 1/kT, where kT is the temperature in energy units. A nontrivial thermo- 
dynamic interpretation is only possible for the nonequal dimer energies (weights), 
Ei ^ E2 (ti 7^ ti). If ti — t2, the energy of any closed-packed dimer configuration is 
the same and there is no thermodynamics. The thermodynamic functions can all be 
deduced from the free energy. There is no phase transition in the dimer model on 
a rectangular lattice [2,3,6,8]. In particular, the specific heat is a smooth function 
of temperature getting its maximum at (kT / AE) max = 0.477572, at which point 
(C/k) max = 0.169375 [8], where AE = \Ei — E 2 \, C/k is the dimensionless specific 
heat, k is the Boltzmann constant. 

Let us now consider the closed-packed dimer model on the standard hexagonal 
lattice. There are three different weights, £1, i2, ^3, connected to each site. The 
hexagonal lattice is equivalent to the so-called brick lattice which can be obtained 
from a rectangular lattice net by the elimination of some of the vertical bonds [6,7]. 
Therefore, the fermionic integral for Q for the hexagonal lattice can be deduced 
from (4.3) and/or (4.4). Another possibility, which we follow below, is to interpret 
the hexagonal lattice (yet on a rectangular net) in a diagonally layered fashion, as 
it was done for the 2D Ising model in [12]. Within such interpretation, omitting 
further details, we obtain the fermionic integral for the partition function of the 2D 
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hexagonal dimer model in a particularly simple form [8]: 

Q = II dc *mndc mn exp ]T [c mn (t! C mn - t 2 C* m _ ln - t 3 Cn-l)] > ( 5 ' 8 ) 
mn mn 

where c mn , c^ n are anticommuting Grassmann variables. Assuming the periodic clos- 
ing conditions for fermions, which in the given case is the boundary approximation, 
we then pass to the momentum space by the simplest Fourier substitution: 

r - 1 V r p i^-(mp+nq) * _ 1 * -j^-(mp+nq) /r q\ 

mn ~ m Z^/ pq ' rnn ~ ry^ pq ' \°-^J 

V Li pq V Li pq 

The integral (5.8) becomes: 

q = /n^^Etv;^-^ 4 ^-^^')] . ( 5 - 10 ) 

pq pq 

with evident explicit solution for Q. Respectively, the free energy per site for the 
hexagonal lattice appears in the form: 

-0/hex=2/ j jUi^[^- t ^ P - t ^ iq ]- ( 5 - U ) 

0^ ' 

In the thermodynamic interpretation, the weights are: t a = exp(e a ), where e a = 
f3E a are the dimensionless energies, the true energies of dimers are — E a , and (5 = 
1/kT. A nontrivial thermodynamics is only possible if the dimer energies (weights) 
are not all equal to each other. From (5.11) it follows that there is the phase 
transition, with frozen low-temperature phase, in the dimer model on the hexagonal 
lattice. To be definite, let t\ < t 2 < t 3 (or —E 3 < —E 2 < —Ex), the critical 
temperature is then given by the condition ti + £2 — ^3 = 0. At low temperatures, 
ti+t 2 < £3, the system is frozen. In this phase (T < T c ) all dimers are captured at the 
bonds with the highest weight (the lowest energy) and the specific heat identically 
equals zero. In the high-temperature phase, t\ + t 2 > t$, the system is unfrozen. 
In this phase (T > T c ) the specific heat exhibits the |T — T c |~5 singularity as T 
approaches T c from above. In the T > T c phase the specific heat is [8]: 

1 A(e 1 ,e 2 ,e 3 ) 



C 
J 



HEX 71 U h +h+ t3 ) + h+ t3 ) ( tl _ h + f 3 ) ( tl + t 2 -t 3 ) 



(5.12) 



where the numerator is a nonsingular function of temperature: A (e\, e 2 , £3) = (e\ — 
£2)(£i — £3) t\ + {e 2 — Si){e 2 — £3) t 2 , + {e 3 — ei){e 3 — e 2 ) t\ . The unusual phase transition 
in the hexagonal dimer model with frozen T < T c phase was first mentioned by 
Kasteleyn in 1963, for a recent comprehensive discussion see [6]. The anisotropic 
specific heat in the form (5.12) is cited here from [8]. 
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6 Conclusions 



We have discussed some aspects of the fermionic interpretation of the 2D closed- 
packed dimer problem. In its original state, the 2D dimer model is rather a discrete 
combinatorial problem. Introducing the fermionic variables into the original parti- 
tion function, we add, effectively, some amount of new degrees of freedom to the 
original discrete sum. In a sense, this extension of the inherent symmetry just makes 
the dimer problem tractable analytically. The dimer partition function can be trans- 
formed into a fermionic Gaussian integral even for the inhomogeneous distribution 
of dimer activities over the lattice bonds. For simple homogeneous lattices, the 
analytic solutions then follow by passing to the momentum space for fermions. 
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